EXAFS, XRD and RMC studies of an Amorphous GasoSeso Alloy Produced by 

Mechanical Alloying 
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The local atomic order of an amorphous GasoSeso alloy produced by Mechanical Alloying (MA) 
was studied by the Extended X-ray Absorption Fine Structure (EXAFS) and X-ray Diffraction 
(XRD) techniques and by Reverse Monte Carlo (RMC) simulations of its total x-ray structure 
factor. The coordination numbers and interatomic distances for the first neighbors were determined 
by means of EXAFS analysis and RMC simulations. The RMC simulations also furnished the partial 
pair distribution functions G^_G^{r), Gg^so('') and Gs^§,{r). The results obtained indicated that 
there are important differences among the local structure of the amorphous GasoSeso alloy produced 
by MA and those of the corresponding crystals, since there are Se-Se pairs in the first coordination 
shell of the amorphous alloy that are forbidden in the GasoSeso crystals. 

PACS numbers: 61.10.Ht, 61.10.Eq, 61.43. Bn, 61.43.Dq, 05.10.Ln, 81.05.Gc 
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In the recent years there has been an increase in the 
number of applications related to nonlinear optical ma- 
terials. However, the desired properties concerning this 
kind of applications, such as optical homogeneity, laser 
damage threshold, stability of the compound upon expo- 
sure to laser beam, ease of fabrication, improved mechan- 
ical strength and the possibility of making large crystals 
are difficult to find in a single material. Some of the ma- 
terials that can be used for nonlinear applications, like 
silver gallium selenides, zinc germanium phosphides and 
thallium arsenic selenides, do not fulfill all of these re- 
quirements, limiting severely their efficiency and applica- 
bility. Thus, there is a high necessity of developing new 
materials with a higher level of performance and more 
cost effective characteristics. Gallium selenide (GaSe) 
has a number of interesting properties for electrical and 
nonlinear optics applications. It transmits in the wave- 
length range varying from 0.65 to 18 /xm and its optical 
absorption coefficient remains below 1 cm~^ throughout 
the transparency range. It has the possibility of con- 
verting sum and difference frequencies Q, 0] ■ It has also 
been reported to be used in making a number of devices 
like MOSFET, IR detector. Solar Cell, compound semi- 
conductor, etc. in crystalline form while in amorphous 
form, it is a potential candidate for optical memory type 
applications P, Q]- Crystalline GaSe is a semiconduc- 
tor of the III- VI family like GaS and InSe and it has 
a layered graphite type structure with a fourfold layer 
in the sequence Se-Ga-Ga-Se. The crystal cleaves very 
easily along the layers 0. At room temperature, the 
layers are bound by weak van der Waals-type interac- 
tions. The weakness of this interaction explains the ex- 
istence of a number of polytypes 0. GasoSeso alloys 
can be prepared by the melting, vapor deposition and 
molecular beam epitaxy techniques H, S 0- These 
techniques have had very limited success because they 
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FIG. 1: Fourier transformation of experimental EXAFS spec- 
tra: a) at the Ga K-edge and b) at the Se K-edge. 



do not have control over the kinetics and morphology. In 
addition, due to the low melting points of the elemental 
Ga (30°C) and Se (217°C) and the high vapor pressure 
of Se above 600° C it is difficult to obtain Ga-Se alloys at 
specific compositions. On the other hand, the mechani- 
cal alloying (MA) technique |^ can be used to overcome 
these difficulties since the temperatures reached in MA 
are very low, what reduces reaction kinetics, allowing the 
production of poorly crystallized or amorphous materials 
jlfl. Illl I12L even if the constituents of the alloy have 
low melting points, as it is in the case of gallium and 
selenium. 
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Amorphous GasgSeso (a-GasoSeso) was prepared by 
MA starting from high purity elemental powder selenium 
(Alfa Aesar, 99.9999% purity, 150 ^m) and scraped in- 
gots of gallium (Aldrich, 99.999% purity) with nominal 
composition GasoSesQ. The mixture was sealed together 
with several steel balls into a cylindrical steel vial under 
argon atmosphere. The ball-to-powder weight ratio was 
10:1. A Spex Mixer/Mill model 8000 was used to perform 
MA at room temperature. The mixture was milled for 15 
h. A ventilation system was used to keep the vial tem- 
perature close to room temperature. The composition 
of the alloy was confirmed by an energy dispersive spec- 
troscopy (EDS) measurement and impurity traces were 
not observed. The alloy produced was investigated by 
extended x-ray absorption fine structure (EXAFS) and 
x-ray diffraction (XRD) techniques and also by reverse 
Monte Carlo simulation (RMC). The EXAFS measure- 
ments were carried out on the D04B beam line of LNLS 
(Campinas, Brazil), using a channel cut monochromator 
(Si 111), two ionization chambers as detectors and a 1 
mm entrance slit. All data were taken at room tempera- 
ture in the transmission mode. The EXAFS oscillations 
x{k) at both K edges, after the standard data reduction 
procedures using Winxas97 software Q, were Fourier 
transformed (FT) using a Banning weighting function 
within the ranges 3.8 - 14.3 A-i for the Ga, and 3.4 - 
14.3 A^^ for the Se edge. They can be seen in Fig.^ Raw 
spectra were filtered by Fourier transforming fc'^x(^) (Ga 
edge) and kx{k) (Se edge) into r-space (Fig. 01 and trans- 
forming back the first coordination shells (1.20 - 2.85 A 
for the Ga edge and 1.13 - 3.0 A for the Se edge). Filtered 
spectra were then fit by using Gaussian distributions to 
represent the homopolar and heteropolar bonds . The 
amplitude and phase shifts relative to the homopolar and 
heteropolar bonds needed to fit them were obtained from 
ab initio calculations using the spherical waves method 
and by the FEFF software. Figure [21 shows the ex- 
perimental and the fitting results for the Fourier-filtered 
first shells on the Ga and Se edges. Structural param- 
eters extracted from the fits are listed in Table It is 
interesting to note that the very good fits shown in Fig. [3 
were achieved only when Se-Se pairs were considered in 
the first shell. This fact indicates that the local struc- 
ture of a-CasoScso produced by MA is different from its 
crystalline counterparts as none of the known stable crys- 
talline GasoSeso structures contains Se-Se bonds. 

The XRD measurements were carried out at the BW5 
beamline |0| at HASYLAB. All data were taken at room 
temperature. The energy of the incident beam was 121.3 
keV (A = 0.102 A). The structure factor 5(i^) (Fig.El full 
line) was computed from the normalized intensity I{K) 
according to Faber and Ziman and it was modeled 
by reverse Monte Carlo simulations. This technique is 
described in details elsewhere and it has 

been used as a method for structural modeling based di- 
rectly on experimental data. There are several papers 
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FIG. 2: Fourier-filtered first shell (full line) and its simulation 
(squares) for a-GasoSeso at the (a) Ga K edge, (b) Se K edge. 




FIG. 3: Experimental (full line) and simulated (squares) total 
structure factor for a-GasoSeso. 



[3 m in il m Hi reporting structural studies of 
amorphous alloys by RMC. Simulations were carried out 
by the RMC program available on the Internet . Cu- 
bic cells contained 1600 and 12800 atoms, and the average 
density was po = 0.03907 ± 0.0005 atoms/A^. This value 
was found from the slope of the straight line {—Airpor) 
fitting the initial part (until the first minimum) of the to- 
tal G(r) function |27|. The minimum distance of atoms 
was also extracted from G{r) and fixed at 2.18 A. All the 
simulations were performed considering atoms randomly 
placed in the cubic cells as starting configurations. Then 
the following series of simulations were carried out: 
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FIG. 4: Gg^g,(r), G%^_g{r) and Gf^^iiir) functions ob- 
tained from the RMC simulations. 

1. Hard sphere simulation without experimental data 
to avoid possible memory effects of the initial con- 
figurations in the results. 

2. 'Unconstrained' runs (i.e. when experimental data 
were 'switched on'). These runs led to three es- 
sentially identical partial pair correlation functions 
and partial structure factors which can be consid- 
ered as linear combinations of the 'true' partial 
quantities. It is to be mentioned that as neither the 
size nor other a priori information can distinguish 
between Ga and Se atoms at this step no adequate 
coordination numbers can be obtained. 

3. 'Constrained' runs. The experimental S{K) was 
fit by using EXAFS coordination number values as 
constraints. Comparison of experimental (full line) 
and calculated (squares) structure factors for the 
latter case is shown in Fig. O and the partial pair 
correlation functions (pcf 's) are given in Fig. 01 

Finally the whole series of calculations was repeated 
from the very beginning with the difference that dur- 
ing the 'constrained' run random steps resulting in non- 
zero Se-Se first coordination number were rejected. It 
is important to note that if Se-Se pairs are forbidden as 
first neighbors simulations did not converge, reinforcing 
the results obtained by EXAFS analysis. The position 
of the first and second peak are 2.42 A and 3.89 A in 
all of the pcf's corresponding to a mean bond angle of 
107° for the four bond types (Ga-Ga-Ga, Se-Se-Se, Ga- 
Se-Se, Ga-Ga-Se) that can be directly derived from the 
pcf peak positions. As this is very close to the value 



TABLE I: Structural parameters obtained for a-GasoSeso. 



EXAFS 





Ga K 


-edge 


Se K-e 


dge 


Bond Type 


Ga-Ga 


Ga-Se 


Se-Ga 


Se-Se 


N 


1.3 


2.4 


2.4 


1.3 


r (A) 


2.38 


2.45 


2.45 


2.37 


(AxlO-2) 


1.58 


0.545 


0.545 


1.77 


RMC 


Bond Type 


Ga-Ga 


Ga-Se 


Se-Ga 


Se-Se 


N 


1.2 


2.5 


2.5 


1.3 


r (A) 


2.42 


2.42 


2.42 


2.42 




GasoSeso 


compound " 






Bond Type 


Ga-Ga 


Ga-Se 


Se-Ga 


Se-Se 


N 


1 


3 


3 


6 


r (A) 


2.44 


2.45 


2.45 


3.75 




GasoSeso 


compound ' 






Bond Type 


Ga-Ga 


Ga-Se 


Se-Ga 


Se-Se 


N 


1 


3 


3 


6 " 


r (A) 


2.39 


2.47 


2.47 


3.74 



"Space group P63/MMC. 
''Space group P-6M2 

"^The trigonal crystal of space group R3M has 4 Se-Se pairs. 



describing perfect tetrahedral coordination (109.5°) and 
A^Ga-Ga + A^Ga-Sc and A^sc-Sc + ^Sc-Ga are both close to 
4 it is evident to assume that ball milled a-GasoSeso has 
a tetrahedral structure with a definite tendency to form 
homopolar bonds. 

The difference of Ga-Ga and Ga-Se bond lengths in 
the crystalline modifications is not greater than about 
0.08 A (see Table and they are also quite close to 
the value of 2.35 A found recently for a-Se [l^. As the 
spatial resolution of diffraction experiments is equal to 
Tr/ATinax the S{K) factor should be measured at least up 
to 40 A~^ or further to get more detailed information 
on the first coordination shell. It should also be men- 
tioned that due to the value of neutron scattering lenghts 
(&Sc = 7.970 fm, bca = 7.288 fm) neutron diffraction data 
would give essentially the same information. Other tech- 
niques used to obtain information at the level of pcf's are 
either prohibitively expensive (neutron diffraction with 
isotopic substitution) or yield limited spatial resolution 
due to the low ATmax value available (anomalous X-ray 
scattering) . 

In summary the local structure of ball milled a- 
GasoSeso was investigated experimentally with EXAFS 
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and high energy X-ray diffraction. EXAFS analysis led 
to the following conclusions: the average first coordina- 
tion number in a-GasoSeso is close to 4, Ga and Se local 
environments are similar and Se-Se bonding is significant. 
All of these findings were checked and confirmed by RMC 
study of diffraction data: it was possible to obtain a good 
fit with coordination constraints close to the EXAFS val- 
ues while runs without Se-Se first neighbors led to a bad 
agreement between model and experiment. Mean bond 
angle calculated from diffraction data is 107° indicating 
that a-Ga5oSe5o has a tetrahedral local structure. 

The present study illustrates how complementary in- 
formation obtained by different experimental techniques 
can be combined within the frame of reverse Monte Carlo 
simulation. We believe that this is a useful and efficient 
way of modelling disordered materials especially in cases 
when traditional methods (e.g. neutron diffraction with 
isotopic substitution) are not available. 
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